Spin dynamics in a dissipative environment: from quantal to classical 
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We address the problem of spin dynamics in the presence of a thermal bath, by solving exactly 
the appropriate quantum master equations with continued-fraction methods. The crossover region 
between the quantum and classical domains is studied by increasing the spin value S, and the 
asymptote for the classical absorption spectra is eventually recovered. Along with the recognized 
relevance of the coupling strength, we show the critical role played by the structure of the system- 
environment interaction in the emergence of classical phenomenology. 



PACS numbers: 03.65.Yz, 05.40.-a, 75.50.Xx, 76.20.+q 

Quantum mechanics is one of the most subtle and pow- 
erful theoretical constructions of the human mind. Un- 
derstanding its implications, relation with other theories, 
and domain of validity has captivated scientists since its 
advent. This domain has been slowly expanding from the 
traditional of atoms and molecules, to condensed matter 
systems (solids and liquids), and more recently by studies 
of decoherence, quantum analogs of classical effects (e.g., 
chaos), and the quantum-to-classical transition 0. 

These studies have also brought an increasing aware- 
ness of the role of the environment. Thus, the field of 
open quantum systems deals with systems consisting of a 
few relevant degrees of freedom coupled to the surround- 
ing medium, which has a large number of constituents 
(photons, phonons, electrons, nuclei, etc.). The coupling 
produces dissipation, fluctuations, and decoherence; it 
also enables the system to interchange energy and corre- 
lations with the bath and relax to equilibrium [2| ■ Besides 
its basic interest, the above generic conditions make this 
topic relevant in various areas of physics and chemistry. 

Spins constitute one of the most paradigmatic quan- 
tum systems due to their discrete and finite energy spec- 
trum. Their dynamics is also special and rich because of 
the underlying commutation relations [Si , Sj] = ie^fciSfc. 
Naturally, it is important to take into account environ- 
mental effects in spin problems, and this has led to sev- 
eral theories of spin relaxation. To deal rigorously with 
quantum dissipative systems, however, is a difficult task. 
Path integral propagators, quantum Langevin or master 
equations can typically be solved in a few simple cases: 
free particle (or in a uniform field), harmonic oscillator 
2], two-state systems (e.g., S = 1/2 spins) Q, etc. 

The continued-fraction method, devised originally for 
classical Brownian motion has been successfully 
adapted to solve master equations for some quantum sys- 
tems 0, 0, 0, 0] . Here we shall apply this technique to 
a spin with arbitrary S weakly coupled to a dissipative 
bath, and monitor its intrinsic dynamics via spin reso- 
nance. We investigate how the approach to the classical 
results takes place (out of reach of previous exact meth- 
ods due to their limitations in S). We focus on the effects 
of the environment, not only of the coupling strength, 



but also of the structure of the spin-bath Hamiltonian. 
Usual studies of open quantum systems overlook the lat- 
ter and adopt the simplest bilinear interaction. We con- 
sider two models with a solid-state motivation: coupling 
to electron-hole excitations, actually linear in S, and to 
phonons, an even polynomial in S. We find that the ap- 
proach to the classical results depends qualitatively on 
the coupling structure. This is specially critical for the 
uniformity of the convergence in the different frequency 
sectors of the spin absorption spectra. The problem is 
not merely academic; large-spin molecular clusters are in 
the focus, while magnetic nanoparticles provide a natural 
classical limit 0. Thus, our results could also help dis- 
criminating different proposed couplings in those systems 
and hence ascertain the microscopic origin of dissipation. 

Let us start with the Hamiltonian of a spin S coupled 
to a bosonic bath (linearly in the bath variables) 

Wtot = H(S) + E q V q F q (S) (o+ + a _,) + H h . (1) 

Here Ti-(S) and Hh = J2 q LU q a q a q are the spin and 
bath Hamiltonians, F q (S) the spin-dependent part of 
the interaction, and V q coupling constants. For sys- 
tems with discrete and finite spectrum, it is very con- 
venient to introduce the Hubbard (level-shift) operators 
X™ = \n)(m\. Any operator F can then be expanded 
in this basis F = Y, nm F nm X™, with F nm = (n\F\m). 
When coupled to the environment, the spin is not in a 
pure state and it needs to be described by its density ma- 
trix g. Its matrix elements are given in this framework 
by g mn = (X?). 

Many problems in quantum optics, magnetism, or 
chemical physics involve weak system-bath coupling 0|. 
Then, the dynamical equation for g can be obtained by 
perturbation theory. In the Hubbard formalism, using 
the S z eigenstates as basis, S z \m) — m\m), one finds for 
H = Hd (S z ) — B ■ S the density-matrix equation 0, 0] 

= iA nm X™ + (i/2) J B + (^ m A™+ 1 -4 1 -iX™ 1 ) 
+ (i/2)B_ (C-i X™- 1 - 4 AT +1 ) + K • (2) 

The A„ m = e n — e m are the frequencies associated with 
the m — > n transition, being e m the levels of the diagonal 
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part of the spin Hamiltonian (including B z ). The circular 
components of the transverse field are B± = B x ± \B y 
and the £ m = [S(S+ 1) — m(m + 1)] 1 / 2 are ladder factors. 

The first three terms in Eq. give the unitary evolu- 
tion of the isolated spin in the Heisenberg representation. 
The relaxation term i?™ incorporates the effects of the 
bath and has a non-Markovian (history dependent) form 

K = - f dr { K(T-t)F( T )[F,XZ} 

-K(t-T)[F ) X?]F(T)}. (3) 

Here the operators without time argument are evaluated 
at t whereas F(r) = £ n , m , F n > m ,X%j (r). The memory 
kernel K. is given in terms of the spectral density of bath 
modes, J(w) = \ J2 q l^gl 2 — an d bosonic occu- 
pation numbers, — (e"/ T — by 



JC{t) = / — J( w )[n„e + " 

JO ^ 



+ (n w + l)e- iwT ] . (4) 



To second order in the interaction, and not too strong 
transverse field, the retarded time dependences X™ (r) 
can be determined by the dominant term in the conser- 
vative evolution X™(r) ~ e~ iA — X™(t). Inserting 
such X™ (r) in i?™, only operators evaluated at t re- 
main and non-Markovian features effectively disappear. 
Then, the coefficients of the X™, include, along with 
the coupling matrix elements F nm , the relaxation rates 
W n \ m = W(A nm ), with the universal rate function asso- 
ciated to the kernel W{A) = Re^dr e - iAr JC(r)]. 

We shall consider in the sequel the following family of 
couplings F: linear in S± = S x ± iS y but allowing for 
^-dependent "coefficients" v(S z ): 



F(S)= V+ [v(S z ), S-]++ V -[v{S t ), 5 



+J + 



(5) 



Here rj± are some scalars ensuring F + — F while 
[A, B] + = A B + B A. Then, the matrix elements F nm — 



{n\F\m) read F nm = L n 



-1 $n,r, 



where L m . m > = r] + [v(m) + v(m')]£ m , m > an d £ m , m ±i = 
[5(5 + I) — m(m ± f)] 1 ^ 2 - The relaxation term for 
these couplings acquires the following (Redfield) form 
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\L m +l,m\ W m _|_i| T 



J rn,m — 1 
f W, 



m|m+l) -^"^+1 • (6) 



Inserting this i?™ in Eq. J2J) we get the master equa- 
tion for our problem within a fully quantum treatment 
(no phenomenological relaxation is introduced, no pre- 
conceived form of the equation is assumed). Recall fi- 
nally that handling the spin precession requires to solve 



the full density-matrix equation, because it involves off- 
diagonal elements, and it is not captured by a (Pauli) 
master equation for the level populations (g m m)- 

We mentioned the difficulties to solve our models 
for quantum dissipation and that the continued-fraction 
method (a relative of the recursion method and Lanc- 
zos tri-diagonalization) has been applied to several quan- 
tum systems 0, S 0- For spin problems ^3] one 
starts writing the master equation compactly as X™ = 

J2n'm' Qn,rf X n' witn Th' = Th — 1, U, Tl+1 and to' = 

m — 1, m, m + 1. To convert this 2-index differential 
recurrence into a I-index one, we introduce appropriate 
(25 + l)-vectors, c„, and (25 + I) x (25 + l)-matrices, 
Qn.n'j with components and elements |12( 



(7) 



obtaining c n = Q„, n _i c„_i + Q„, n c» + <Q> n ,„+i c n+:L . 
In this I-index form the recurrence can be tackled by 
(matrix) continued-fraction methods Q yielding the so- 
lution of the master equation (0) . We have then the full 
density matrix g mn — (X™) — (c n ) m and any observ- 
able (magnetization, susceptibilities, etc.) can be com- 
puted from the trace formula (A) = Tr(gA). The matrix 
associated to the original system, X = Q X , had di- 
mensions (25 + l) 2 x (25 + l) 2 dimculting the handling- 
even of moderate spins (5 < 10). The continued-fraction 
approach replaces it by 25 + 1 problems with matrices 
(25+1) x (25+ 1), allowing to gain some orders of mag- 
nitude in 5 and pursue the classical limit a way longer. 

We now apply the above formalism to the problem of 
spin dynamics in a magnetic-anisotropy plus Zeeman po- 
tential Ti = —D S z —B-S. This Hamiltonian may also be 
viewed as the minimal model for superparamagnets 0|. 
The anisotropy term has two minima at S z = ±5 with a 
barrier at S z =0. The coupling to the environment pro- 
vokes quantum Brownian motion of the spin, which may 
overcome the potential barriers. We consider two basic 
solid-state mechanisms 0, Q : (i) Coupling to electron- 
hole excitations near the Fermi surface (a bosonizable 
bath); then F{S) = |(ry+5_ + 77-5+) [i.e., v(S z ) = const 
in Eq. (JSJ)] while the bath is Ohmic, J {to) = \ u>. (ii) Cou- 
pling to phonons; now v(S z ) oc S z and the environment 
is super-Ohmic, J(oj) — Xlo 3 (in 3D). The rates W n n m / 
required in R™ can be obtained from the spectral density 
by W(A) = J(A) n A + J(-A)(ra_ A + 1), understanding 
J(u> < 0) = 0. We start with the super-Ohmic nonlin- 
ear case, which has received less attention in the context 
of quantum dissipative systems than the Ohmic bilinear 
coupling; we will see that it also has a rich physics. 

The Zeeman term and F oc S± have non-zero matrix 
elements between the states |to), producing transitions 
between them. In an oscillating field, they result in peaks 
in the imaginary part x"(^) of the dynamical susceptibil- 
ity (absorption line-shape) located at the transition fre- 
quencies A mjm+ i = e m -£ m+ i = D(2m+1)+B z (Fig.[TJ|. 
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FIG. 1: Absorption line-shape x"(^) f° r a s Pi n S — 10 with 
D = 0.5 at B z = 0. Thick line: a = D S 2 /T = 5 and spin- 
phonon coupling A = 3 • 10 -8 . Thin lines: effects of halving 
the damping at the same T (solid), and halving T keeping 
A (dashed). Vertical lines: loci of the transition frequencies 
A m ,m+i = D(2m+1). Inset: classical dampingless asymptote 
for a = 1, 2 and 3 (lines), and exact Fokker-Planck results 
for finite Landau-Lifshitz damping (All = 0.003; symbols). 



We thus see that the classical phenomenology looks 
quite different from the multi-peaked structure of the 
quantum case. This poses the following questions: 
(i) How does quantum mechanics manage to join those 
two behaviors? and (ii) which are the main factors de- 
termining the way in which the classical phenomenology 
emerges? We now try to answer these questions by solv- 
ing the density-matrix equation (J2J for increasing 5* and 
getting as close as possible to the classical domain. 

Recall, however, that limiting procedures in physical 
problems (e.g., lattice to continuous limit in field theo- 
ries, thermodynamic limit in statistical mechanics, etc.) 
require to define clearly: (i) which quantities are kept 
constant when taking the limit and (ii) which scaled vari- 
ables are needed to monitor the evolution. We fix the 
reduced anisotropy and field parameters a = DS 2 /T 
and £ = SB/T. At constant T this entails keeping the 
anisotropy-barrier height and amount of Zeeman energy 
constant (and hence finite) while introducing more lev- 
els with S (the spacing decreases as A ~ 1/S). As for 
the scaled quantities, guided by the classical result JHJ, 
we use x/xo with xo — S(S + 1)/T (corresponding to 



The transitions at the potential wells (\m\ ~ S) corre- 
spond to the largest frequencies (A ~ 2DS at B z = 0), 
while those near the barrier top (m ~ 0) appear at low f2 
(~ D). Then, going from high to low f2, the intensity of 
the peaks decreases, as they involve transitions between 
higher levels, which are thermally less populated. 

The peaks have finite width and height due to the 
damping A and the temperature, as the interaction with 
the bath "blurs" the spin energy levels. Thus, a lowering 
of A or T makes the peaks narrower and higher (phe- 
nomenology akin to that of a damped oscillator). There 
is an extra narrowing of the low Q peaks, because the 
spin-phonon coupling F ~ S Z S± leads to an effective 
damping decreasing with m [A e g ~ A(2m± 1) ]. This en- 
ters in i?™ via the modified ladder factors |L m ,m±i| 2 ~ 
(2m ± 1) 2 ^„ m ±i and it is the spin analogue of position- 
dependent damping in translational Brownian motion. 

Next, let us briefly discuss the corresponding classical 
behavior. The actual line-shape will depend on the phe- 
nomenological relaxation model considered (Bloch equa- 
tions, Landau-Lifshitz, etc). Nevertheless, the result in 
the limit of zero damping is universal yjj |l4j 
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x"(fi) = n[1 ~ w n -^ ex p[^/ao 2 ] • (8) 

Here Z is the partition function, f7 a the resonance fre- 
quency at the wells, and a the barrier over T. Physically, 
the anisotropy Hd = —D S z leads to S'z-dependent pre- 
cession frequencies and the ensuing spreading of the spec- 
tral line-shape (inset of Fig.GJ. The population of the dif- 
ferent Sz-orbits changes with T, modifying x"(£l). Note 
that this dissipationless limit provides a good description 
for weak enough coupling in most of the range. 
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FIG. 2: Spectral line-shape x"(fy for o = 1 at B z = 
(see the text for the scaled quantities). The thick dashed line 
is the classical Eq. ©. Top: 5 = 5, 25, 50, and 100 with 
"constant" coupling X/S = 10 -2 . Bottom: fixed S = 50 with 
X/S = 10~ 2 (as in top), 3 ■ 10" 2 , and 10~\ 
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FIG. 3: Line-shape for S — 10 at a — 1 with (i) the phonon- 
coupling model and (ii) bilinear coupling, both with super- 
Ohmic (hybrid) and Ohmic (electron-hole) spectral densities. 



/i 2 /T) and Q/2DS (which tends to fi/fi a ). Finally, we 
also scale the bare coupling strength A with S. Inspecting 
the density-matrix equation we see that the Hamiltonian 
coefficients go as A ~ 1/S while the relaxation ones de- 
crease as XD 2 L 2 A 2 ~ X/S 2 (we include a D 2 dependence 
arising in the coupling to phonons |llj|). Thus, fixing X/S 
we can study the effects of going to large S while main- 
taining the relative "weights" of the conservative and re- 
laxation terms in the quantum master equation. 

Proceeding in this way, we compute the transverse dy- 
namical response for various S (Fig. For moderate 
spins we clearly recognize the quantum features of Fig.^ 
As S is increased more peaks are introduced into the same 
interval Q/O a . Due to their finite width they start to co- 
alesce and a limit curve progressively emerges. However, 
the approach is far from uniform in CI. At low frequen- 
cies the peaks merge slowly with S; they are sharp and 
narrow due to the m-dependent damping associated to 
S z in F ~ S Z S±. This is less relevant at high frequen- 
cies [large m|, SX c s/X c s ~ 4 /(2m ± 1)] and a smooth 
peakless line-shape arises there. For a fixed S, in addi- 
tion, one would expect that larger spin-bath coupling will 
"accelerate" the classical convergence. Figure actually 
shows that the wildly peaked part is then pushed further 
into the low 17 sector and that the "oscillations" around 
the limit curve are reduced. It is remarkable that this 
limit curve is indeed Gekht's classical prediction |JSJ. 

Finite width of the absorption peaks has been essen- 
tial to reconstruct the classical curve. Here it has been 
provided by the coupling to the environment; in other sit- 
uations different broadening mechanisms may contribute 
(lif. The form of the interaction, on the other hand, 
has led to a highly non-uniform approach to the classical 
asymptote. This shows that not only the strength, but 
also the structure of the coupling Hamiltonian can play 
an important role in the approach to the classical regime. 



This can be further supported by comparing with the 
"electron- hole" coupling model, where F ~ S±. To as- 
sess the different contributions we proceed in two steps, 
adjusting A to get the same intensity for the ground-state 
transitions; Fig. [21 First, we go from the phonon-coupling 
F ~ S Z S± to a hybrid model with F ~ S±, but still 
super-Ohmic spectral density. This F greatly tames the 
low-frequency sharp peaks, but some non-uniformity still 
remains, due to W m \ m+ i ~ J (A) ua ~ ^-rnm+l in the re- 
laxation term. Second, we add the Ohmic bath J cx lo to 
the bilinear coupling. Then W m | m+1 ~ A-mm+i ~ const 
and the approach to the classical behavior becomes quite 
uniform in most of the 17 range, in spite of the moderate 
spin value considered (S = 10). 

In summary, we have addressed the problem of spin 
dynamics in a dissipative thermal bath. Solving exactly 
the quantum master equation by a continued-fraction 
method for increasing S has allowed us to approach the 
classical prediction for the absorption spectra. We have 
investigated the effects of the spin-bath interaction on the 
quantum-to-classical crossover. The coupling strength, 
as usual in quantum dissipative systems, accentuates the 
attainment of the classical phenomenology. However, the 
approach is qualitatively affected by the structure of the 
interaction, as illustrated with the different convergences 
in the different sectors of the absorption spectra for two 
important solid-state mechanisms. Although the rele- 
vance of dissipation, specially in mesoscopic systems, is 
amply recognized, only studies of decoherence and ap- 
proach to equilibrium had paid due attention to the struc- 
ture of the coupling Hamiltonian. Here we have shown 
its relevance also in the features of the quantum-classical 
border and in the emergence of classical behavior. 
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